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ABSTRACT 

V Sagittae shows quasi-periodic optical high (soft X-ray off) and low (soft X-ray on) states with the 
total period of ^ 300 days. A binary model is presented to explain orbital light curves both for the 
high and low states as well as the transition mechanism between them. The binary model consists of 
a white dwarf (WD), a disk around the WD, and a lobe-filling main-sequence (MS) companion. In the 
optical low state, the mass transfer rate to the WD is small and the size of the disk is as small as its 
Roche lobe size. In the optical high state, the mass transfer rate to the WD exceeds the critical rate of 
~ 1 X 10~^Mq yr~^, and the WD blows an optically thick, massive wind. Surface layers of the disk are 
blown in the wind and the disk surface extends to the companion or over. As a result, optical luminosity 
of the disk increases by a magnitude because of its large irradiation effect. The massive wind completely 
obscures soft X-rays. This corresponds to the optical high/soft X-ray off state. The transition between 
optical high and low states is driven by an attenuation of the mass transfer from the secondary. During 
the optical high state, the wind from the WD hits the companion and strips off its surface layer. The 
mass transfer from the companion gradually reduces and stops. As the mass supply stops, the WD 
wind weakens and eventually stops. The disk shrinks to a Roche lobe size and the optical magnitude 
drops. This phase corresponds to the optical low/soft X-ray on state. Then a rapid mass transfer 
resumes. The mass of the WD envelope increases and it blows a wind again. This cycle is repeated like 
a limit cycle. During this intermittent wind phase, the WD can grow in mass at the critical rate and 
eventually reach the Chandrasekhar mass limit. This process is called "accretion wind evolution," which 
is a key evolutionary process in a recently developed evolutionary scenario of Type la supernovae. This 
evolutionary process was first confirmed in the LMC supersoft X-ray source RX J0513.9— 6951, although 
it commonly occurs in the supersoft X-ray sources when the mass transfer rate exceeds the critical rate. 
Thus, V Sge is the second example of accretion wind evolution. 



Subject headings: binaries: close — novae, cataclysmic variables 
stars: winds, outflows — X-rays: stars 



stars: individual (V Sagittae) 



1. INTRODUCTION 

V Sagittae is an eclipsing binary system exhibiting op- 
tical high and low states (Herbig et al. 1965). V Sge has 
been intensively observed in various bands from radio to 
X-rays and many groups challenged to interpret its binary 
nature, but this peculiar system still defies explanation 
(e.g., Smak et al. 2001, for recent review and criticism on 
the binary nature). 

There have been many debates on the nature of V Sge. 
Williams et al. (1986) suggested a cataclysmic variable 
(CV) nature from their observations. They claimed that 
the observed changes in the line profiles through the eclipse 
were similar to the rotational disturbance expected from 
the eclipse of an accretion disk. However, Wood & Lockley 
(2000) argued that the variability of the emission lines seen 
by Williams et al. (1986) continue all the way around the 
orbit and is actually unrelated to the eclipse. Robertson, 
Honeycutt, & Pier (1997) also found that no evidence of a 
rotational disturbance in their high-state spectroscopy of 



Ha. 

V Sge has been suggested to be a supersoft X-ray source 
(SSS) (Diaz & Steiner 1995; Greiner & van Teeseling 
1998; Patterson et al. 1998; Steiner & Diaz 1998). 
This suggestion stems from similarities in the orbital light 
curve shapes between V Sge and the LMC supersoft X- 
ray source CAL 87 and in the long-term soft X-ray behav- 
iors between V Sge and the LMC supersoft X-ray source 
RX J0513.9— 6951. Soft X-rays were detected only during 
the optical low states (Greiner & van Teeseling 1998), al- 
though this is not yet fully established in V Sge because 
of their sparse period coverage in X-rays. 

Wood & Lockley (2000) argued that the light curves of 
V Sge can be well fitted by a much simpler model (a pair 
of two stars, see also Mader & Shafter 1997; Lockley et 
al. 1999; Smak et al. 2001) and that no one has yet mod- 
eled the light curves or explained the emission-line profiles 
wfth the SSS model. Instead, Lockley (1998) and Wood & 
Lockley (2000) proposed a colliding wind model of V Sge 
both in the optical high and low states (see, also, Gies et 
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al. 1998; Lockley & Wood 1998; Wood & Lockley 1997). 
However, it should be noted that no one successfully pro- 
posed a model of the transition between the optical high 
and low states. 

The key observational features necessary to be eluci- 
dated here are: (1) V Sge exhibits long-term transitions 
between optical high (y ~ 11 and ~ 180 days) and low 
(V ^ 12 and ~ 120 days) states with the total durations of 
about 300 days (e.g., Simon & Mattei 1999, for the long- 
term behavior). (2) Very soft X-rays are detected only in 
the long-term optical low state (e.g., Greiner & van Teesel- 
ing 1998). (3) Recent radio observations indicate a wind 
mass loss rate of ^ 1O~^M0 yr~-^ (Lockley, Eyres, & Wood 
1997; Lockley et al. 1999). 

These observational characteristics remind us of the 
LMC supersoft X-ray source RX J0513.9— 6951 (hereafter 
RX J0513). RX J0513 shows very similar characteris- 
tics: (1') quasi-regular transitions between optical high 
{V ~ 16.5 and ~ 120 days) and low {V ~ 17.3 and ~ 40 
days) states (e.g., Alcock et al. 1996; Cowley et al. 2002), 
(2') supersoft X-rays detected only in the optical low states 
(e.g. Reinsch et al. 1996, 2000; Schacidt et al. 1993; 
Southwell et al. 1996), and (3') bipolar jets and winds 
with the velocity of ~ 4000 km s~^ (e.g., Crampton et 
al. 1996; Hutchings et al. 2002; Southwell et al. 1996), 
which is much stronger (faster and denser) in the optical 
high states. Recently, a new transition mechanism be- 
tween the optical high and low states of RX J0513 has 
been proposed by Hachisu & Kato (2003a,b). Here we 
try to reproduce the above three observational features of 
V Sge by the same mechanism as in RX J0513. 

In §2, we introduce our binary model for V Sge. The 
orbital light curve modulations in the optical high and low 
states are reproduced in §3. We discuss various arguments 
for/against the WD model in §4. A limit cycle model for 
the long-term light curve behavior is formulated in §5 and 
we present our numerical results in §6. Discussion and 
conclusions follow in §7 and §8, respectively. 

2. THE BINARY MODEL OF V SAGITTAE 

The orbital light curve modulations of V Sge are quite 
different both in shape and brightness between the opti- 
cal high and low states (e.g., Herbig et al. 1965; Mader 
& Shatter 1997; Simon et al. 2002). In this section, 
we will show how to construct our binary models that re- 
produce orbital modulations in the optical high and low 
states. We adopt a binary system consisting of a mass- 
accreting white dwarf (WD), a disk around the WD, and 
a lobe-filling, main-sequence (MS) companion (see an il- 
lustration of Fig. 1). 

2.1. Binary parameters 

In this section, we assume the WD mass of Mwd ~ 
1.25 Mq taking the results of §6 in advance, i.e.. 

Ml = Mwd = 1-25 Mq. (1) 

We adopt the ephemeris by Patterson et al. (1998), i.e., 

i(HJD) = 2, 437, 889.916 + 0.5141967 x E 

-9.2 X 10""^;^ (2) 

at eclipse minima (see also, Lockley et al. 1999; Mader & 
Shatter 1997; Simon 1996b; Smak 1995). 



The secondary mass is determined as follows: Her- 
big et al. (1965) estimated each component mass from 
the observed radial velocities. They assumed that the 
weaker O III fluorescent line correctly traces the mo- 
tion of the primary while the stronger O III fluores- 
cent line follows the motion of the secondary, i.e.. Mi = 
(0.74 ± 0.19) Mq/ sin^ i from Ki ~ (320 ± 28) km s'^ and 
M2 = (2.80 ± 0.65) Mq/ sin^ i from K2 - (85 ± 9) km s'^ 
(Smak et al. 2001), where i is the inclination angle of 
the binary. Radial velocities of these two sharp lines are 
oscillating in the same period of the eclipsing binary and 
in anti-phase each other. Moreover, gamma velocities of 
these two lines are almost coincident. We do not think 
that the other lines such as H and He II are correct tracers 
of binary components because they are broad and their 
gamma velocities are not consistent with those of O III 
fluorescent lines. 

Here, we reexamine the companion mass. We assume, 
following Herbig et al.'s paper, that the weaker O III flu- 
orescent line correctly traces the motion of the primary, 
i.e., Ki ~ 300 - 340 km s'^ (Herbig et al. 1965), but do 
not think that the stronger O III fluorescent line correctly 
follows the motion of the secondary. The reason is as fol- 
lows: When the companion star is strongly irradiated by 
a luminous (hydrogen shell-burning) WD, its surface tem- 
perature is much higher at the irradiated hemisphere than 
that of the other side (see Fig. 1). Probably the stronger 
O III fluorescent line should trace more closely the irra- 
diated hemisphere of the cool component. As a result, 
it is very likely that the observed orbital velocity of the 
stronger O III fluorescent line is a bit slower than the mo- 
tion of the cool component as seen in Figure 1. If this is the 
case, the observed velocity oi K2 = 85 ±9 km s~^ does not 
represent the true radial velocity of the cool companion. 

Figure 2 shows the relation between the orbital radial 
velocities and the MS companion mass derived from Ke- 
pler motion with the period of 0.514 days. This flgure 
shows the secondary mass between Mms = 3.0 and 3.5 Mq 
for Ki 300 — 340 km s~^, and the orbital radial ve- 
locity of K2 = 130 — 110 km s~^. Here, we adopt the 
mass of M2 = AfMs = 3.0 Mq. Then the separation 
is a = 4.375 Rq, the radii of the effective Roche lobes 
are = 1.34 Rq and R* = R2 = 2.0 Rq. The ra- 
dius of 3.0 Mq zero-age main-sequence (ZAMS) star is 
-RzAMS = 2.5 Rq (e.g., Bressan et al. 1993), so that the 
secondary is a main-sequence (MS) star but it is under- 
sized and probably underluminous because of thermal im- 
balance by the Kelvin-Helmholtz timescale mass-transfer. 
We assume TMs,org = 12, 000 K for non-irradiated tem- 
perature of the MS companion. The inclination angle is 
determined by light curve fitting as described below. It 
should be noted here that the thermal timescale mass- 
transfer is realized for the mass ratio M2/M1 > 1 and the 
mass transfer rate can reach ~ 1O~^M0 yr"-*^ or more (e.g. 
Li & van den Heuvel 1997). 

Our binary models for numerical calculation are illus- 
trated in Figure 3. A circular orbit is assumed. We also 
assume that the photospheric surfaces of the WD, the MS 
companion, and the disk emit photons as a blackbody at 
a local temperature that varies with position. The numer- 
ical method adopted here was fully described in Hachisu 
&Kato (2001b). 
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2.2. Flaring-up disk and its irradiation effect 

The irradiation effects both of the disk and of the MS 
companion play an essential role in tlic light curve mod- 
eling. For the basic structure of the accretion disk, we 
assume an axisymmetric structure with the size and thick- 
ness of 



-Rdisk = OiR 



1) 



and 



h = /3i?dii 



isk 



(3) 



(4) 



where i?disk is the outer edge of the accretion disk, is 
the effective radius of the inner critical Roche lobe for the 

WD component, h is the height of the surface from the 
equatorial plane, and zu is the distance from the symmet- 
ric axis. 

Since the orbital light curve modulations of V Sge show 
a nonaxisymmetric shape, we introduce an asymmetric 
configuration of the accretion disk as done by Schandl, 
Meyer-Hofmeistcr, & Meyer (1997). They explained an 
asymmetric feature in the orbital light curves of the LMC 
luminous supersoft X-ray source CAL 87. The cause is 
that the gas stream from the companion hits the edge of 
the accretion disk and makes a vertical spray when the 
mass transfer rate is as high as in the luminous supersoft 
X-ray sources (~ 1 x 10~'^Mq yr~^). We have also intro- 
duced the same type of asymmetry of the accretion disk 
as Schandl et al.'s, i.e.. 
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h \ Clow 



for < ^ < ^1 



for C > fa, 
for C < fa, 



where 



disk 



(5) 
(6) 

(7) 



and Clow is a parameter specifying the degree of asymme- 
try, (f)3 is the starting angle of the vertical spray, at 4>4 the 
vertical spray reaches its maximiim height, this maximum 
height continues until the angle of the elevated rim 
disappear at the angle of (j)2, fa is the starting radius ratio 
from where the vertical spray is prominent, z is the height 
of the disk surface from the equatorial plane for the asym- 
metric case. The azimuthal angle is measured from the 
phase at the WD component in front of the MS compan- 
ion, i.e., from the phase at the secondary eclipse (see Fig. 
1). Here, we assume that the accretion disk is symmetric 
inside zu < 0.8i?disk {fa = 0.8). The other parameters 
of Clow, 4*1 J 4'2, (f'Sj s-iid 04 are determined from the light 
curve fitting. Such an example of the accretion disk is 
shown in Figure 3. Here, we adopt a = 6.0, /3 = 0.24, 
V = 1, Clow = 0.1, (bi = (21/16)^, 02 = (25/16)^, 
03 — (5/16)7r, and 04 = (10/16)7r for configuration (a) 
but a = 1.15, P = 0.4, 1/ = 2, Clow = 0.25, 0i = (23/16)7r, 
02 = (27/16)7r, 03 = (3/16)7r, and 04 = (7/16)7r for con- 
figuration (b). The irradiation efficiency is assumed to be 



Vdk = 0.5, i.e., 50% of the irradiated energy is emitted 
but the residual 50% is converted into thermal energy. 

Recent analysis on the recurrent nova CI Aql 2000 out- 
burst (Hachisu & Kato 2003a) and the LMC supersoft 
X-ray source RX J0513 (Hachisu & Kato 2003b) indicates 
that the irradiation area of the disk is blown in the wind 
and extends up to the companion star or over (a ~ 3 or 
over) as schematically shown in Figure 4. This physically 
means that (1) the surface layer of the accretion disk is 
blown off and accelerated up to. at least, several hundred 
km s~^ like a free stream going outward, and that (2) the 
optically thick part of the disk surface ends somewhere far 
outside the Roche lobe as seen in Figures 3a and 4a. Here, 
wc adopt ly = 1 during the massive wind phase but 1/^2 
in the no massive wind phase (Hachisu & Kato 2003a, b) . 
This tu-square law mimics the effect of flaring-up at the 
rim of the disk by vertical spray as discussed by Schandl 
et al. (1997) for supersoft X-ray sources. 

The accretion luminosity of the disk is also numerically 
included, although its contribution to the F-light is much 
smaller compared with that of the irradiation effect (see 
discussion of Hachisu & Kato 2001a, b; Hachisu, Kato, 
& Schaefcr 2003). The temperature of the outer rim of 
the disk is assumed here to be 10000 K. We have checked 
that the light curves are hardly changed even if 8000 K or 
12000 K is adopted as an outer rim temperature instead 
of 10000 K. 

2.3. Irradiated MS companion 

The MS companion is also irradiated by the WD. As 
already shown by Smak et al. (2001), the irradiation 
effect alone cannot closely reproduce the orbital modula- 
tions. Schandl et al. (1997) adopted a redistribution of 
the irradiated energy with a diffusion angle of 45° to re- 
produce a relatively round shape of orbital modulations 
just before/after the eclipse minima. Here, we also adopt 
the same type of redistributions of irradiated energy, but 
with a bit larger diffusion angle of 60° instead of 45°. The 
irradiation efficiency is the same as that of Schandl et al. 
(1997), i.e., assumed to be tjms = 0.5, i.e., 50% of the 
irradiated energy is emitted but the residual 50% is con- 
verted into thermal energy. The distribution of surface 
temperature is determined by the constraint that the to- 
tal of the emitted energy is equal to the sum of 50% of the 
irradiated energy and the emitted energy with the original 
surface temperature of 12,000 K. 

3. ORBITAL LIGHT CURVES IN THE OPTICAL HIGH AND 
LOW STATES 

Our modeled light curves are shown in Figure 5. To 
reproduce the bright state, we must have a large size of 
the disk, a ~ 6.0, and an expanded WD envelope of 
-RwD.ph ~ 1-7 Rq (TwD,ph ~ 60,000 K), as shown in 
Figure 3a. For the optical faint state, we adopt a much 
smaller photospheric radius of i?wD.ph ^ 0.02 Rq and a 
smaller size of the disk, i.e., a = 1.15, as shown in Fig- 
ures 3b, 6b, and 7b. The 1.25 Mq WD blows winds when 
-RwD,ph ^ 0.07 Rq, so that these two photospheric radii 
are consistent with the optically thick wind theory. 

The photospheric radius of the WD (i?vvD,ph = 1-7 Rq) 
exceeds its critical Roche lobe {Rl = 1.34 Rq) but smaller 
than that of the MS secondary (i?| = R2 = 2.0 Rq). The 
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velocity of the wind is as fast as ~ 1000 km s^^ or more at 
the WD photosphere and much faster than the orbital ve- 
locity (~ 400 km s~^). As a result, the wind is almost free 
stream near the photosphere and hardly affected by the 
rotation effect (or the tidal effect) even if the photosphere 
exceeds the Roche lobe. 

The velocities of winds calculated are as fast as 1000 — 
1500 km s~^ (e.g., Kato k Hachisu 1994) and the wind 
mass loss rate is as large as 10~^ — 10~'^Mq yr~-^ (e.g., 
Hachisu & Kato 2001b). Obscrvationally wind velocities 
are reported to be, for example, 4000 km s~^ during the 
bright state (see, e.g., Koch et al. 1986; Wood & Lockley 
2000). Such a massive and fast wind affects the surface of 
the accretion disk. Because of a large velocity difference 
between the wind and the disk surface, it certainly drives 
the Kelvin-Helmholtz instability at the interface. Only 
the surface layer of the disk is dragged away outward with 
the velocity at least several hundred km s~^ like a free 
stream moving outward. This surface free stream is op- 
tically thick near the original disk but becomes optically 
thin somewhere outside because of geometrical dilution ef- 
fect. We regard the transition place from optically thick 
to thin as the edge of the extended disk (see Fig. 4). In 
this sense, we should understand that the outer rim shown 
in Figure 3a is not the matter boundary but a transition 
from optically thick to thin regions of the tenuous disk 
surface flow. 

It should be noted that a high density part of the disk 
(dark-hatched region in Fig. 4) is hardly changed by this 
relatively tenuous disk surface free stream because the in- 
ternal density of the disk is much denser than that of 
the WD wind. The wind mass loss rate reaches about 
1 X IQ-^Mq yr-i and its velocity is ~ 1000-4000 km s^^, 
so that the density of the wind is estimated from the 
continuity (Mwind = Anr^pv) to be about lO^^^g cm^"^ 
at the distance of 3 Rq from the center of the WD. On 
the other hand, the density of the accretion disk is about 
1 X 10~^g cm~^ at the same radius. Here, we assume the 
standard accretion disk model (Shakura & Sunyaev 1973) 
with a mass accretion rate of ^ 1 x 1O~'''M0 yr^^. 

In the faint state, on the other hand, small WD radii and 
disk sizes are required as illustrated in Figure 3b. As sug- 
gested by Smak et al. (2001), the size of the disk exceeds 
the Roche lobe, that is, a = 1.15, as shown in Figures 3b, 
6b, and 7b. The edge of the disk still contributes to the 
total light even in the eclipse minima (see Fig. 6b) and, 
as a result, the eclipse is not total but partial. 

The orbital light curves both in the high and low states 
vary from night to night (e.g., Herbig et al. 1965; Madcr 
& Shaffer 1997; Lockley et al. 1997). These variations 
are attributed to the dynamical changes of the disk edge 
shape, which are frequently observed in supersoft X-ray 
sources. In this sense, our obtained shapes of the disk 
both in the high and low states are not rigid but dynami- 
cally variable. 

The calculated color variations are roughly consistent 

with the observations, if the color excess is E{B — V) ^ 0.3 
(e.g., Mader & Shaffer 1997; Patterson et al. 1998), i.e., 
for the faintest state, [B — V)o ^ 0.0 against {B — V)c ^ 
—0.3, and, for the brightest state, {B — V)o ~ —0.1 against 
{B — V)c ~ —0.4, where {B — V)c is the calculated color 
(see Fig. 5) and {B — V)o the observed color (taken from. 



e.g., Herbig et al. 1965; Simon et al. 2001; Smak et 
al. 2001). Then, the distance to V Sge is estimated to be 
D = 3.0 kpc with (m-M)y = 13.32 and E{B-V) = 0.30. 

4. VARIOUS ARGUMENTS FOR/ AGAINST THE WHITE 
DWARF MODEL 

The most serious criticisms on the WD model of V Sge 

have been proposed by Smak et al. (2001). They excluded 
the model of a disk around a compact primary component 
mainly because such a model cannot reproduce the orbital 
light curves in the faint phase. However, their conclu- 
sions are derived based on the assumptions that (1) the 
accretion disk is geometrically thin, that (2) their irradi- 
ation effect does not include the diffusion effect between 
the irradiated and non-irradiated regions of the secondary 
surface, and that (3) the size of the disk is limited within 
0.9 times the Roche lobe size. As already shown in the 
previous subsection, we arc able to reproduce the orbital 
light curves if these three assumptions are relaxed in V 
Sge. 

As for assumption (1), Schandl et al. (1997) have shown 
that an elevated edge of the disk reproduce the light curve 
of the supersoft X-ray source CAL 87, which is very simi- 
lar to the faint state of V Sge (e.g., Diaz & Steiner 1995; 
Greiner & van Teeseling 1998; Patterson et al. 1998; 
Simon 1996a). As for assumption (2), Schandl et al. 
(1997) also showed that a diffusion effect on the secondary 
surface is important to explain a wide primary eclipse (or 
a round shape just outside the primary eclipse). The same 
round shape just outside the primary eclipse was also ob- 
served in the recurrent nova CI Aql 2000 outburst and 
such a round shape can be reproduced by the same effect 
(Lederle & Kimeswengcr 2003). 

As for assumption (3), the mass transfer rate becomes 
as large as the Eddington limit for white dwarfs in the 
optical low (faint) state, i.e., 

47rcGMwD RwD , t;,.r 1 /^N 

MEdd = — ^,7 ~ 1 X IO-^Mq yr-i, (8) 

as will be shown in the following section. In such high ac- 
cretion rates, a part of matter overflows the inner critical 
Roche lobe. However, the size of the disk, i.e., a = 1.15, 
is still smaller than the outer critical Roche lobe. 

Hoard, Wallerstein, & Willson (1996) found that the X- 
ray flux from V Sge shows modulation with periodicities 
approximately equal to a half, a third, and a fourth of the 
orbital period but not with the orbital period itself. Smak 
et al. (2001) counted it one of the evidences against the 
model in which the hotter component is a source of soft 
X-rays. In our model, however, the hot component (white 
dwarf) is blocked at least twice, i.e., by the edge of the disk 
and by the companion as seen in Figure 3. Furthermore, 
in some configurations of the disk like in Figure 3, this 
occultation of the hot component by the disk edge contin- 
ues from binary phase ~ 0.4 to 0.7, i.e., ~ 0.3 orbital 
periods, which may explain a modulation frequency of a 
third of the orbital period. 

5. A LIMIT CYCLE MODEL OF V SAGITTAE 

The long-term optical behavior of V Sge shows a quasi- 
regular transition between high and low states (e.g., Simon 
& Mattel 1999, for a recent summary of the long-term 
variations). In this section, we reproduce the long-term 
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optical variations based on an optically thick wind model 
of mass-accreting white dwarfs (Hachisu & Kato 2001b; 
Hachisu et al. 1996, 1999a,b; Kato 1983; Kato & Hachisu 
1994). 

5.1. Winds of mass- accreting WDs 

An essential point of our wind model is that WDs blow 
a fast and massive wind when the mass accretion rate ex- 
ceeds the critical rate, i.e., 

'MwD 



Mace > Mcr « 0.75 X 10"^ 



M, 



© 



-0.4 Moyr-i, (9) 



for X = 0.7 and Z = 0.02 (see, e.g., Hachisu & Kato 
2001b, for other cases of hydrogen content and metal- 
licity). The wind velocity is as fast as or faster than 
I'wind ~ 1000 km s~^ and the wind mass loss rate is ap- 
proximately given by 

Mwind « Mace - M„, (10) 
and usually as massive as ~ 1 x 1O~^M0 yr~^. 

5.2. Mass stripping effects 

When massive winds collide with the surface of the com- 
panion, the surface is shock-heated and ablated in the 
wind. We estimate the shock-heating by assuming that 
the velocity component normal to the surface is dissipated 
by shock and its kinetic energy is converted into the ther- 
mal energy of the surface layer. We regard that gas is 
ablated and lost from L3 point (L3 is the outer Lagrange 
point near the MS companion) when the gas gets the same 
amount of thermal energy as the difference of the Roche 
potentials between the MS surface and L3 point. Then the 
mass stripping rate is given by 



GM 



(?!'L3 - ' 



>MS) 



■ Mstrip = -t^^Mwind • VeS ' Oiq), (H) 



where M = A/wd + Mms; o is the separation of the binary, 
0MS the Roche potential at the MS surface, (j)L3 means 
the Roche potential at L3 point near the MS companion, 
both of which arc normalized by GM/a, rjca the efficiency 
of conversion from kinetic energy to thermal energy by 
shock, g{q) is the geometrical factor of the MS surface hit 
by the wind including the dissipation effect (only the nor- 
mal component of the wind velocity is dissipated), and 
g{q) is only a function of the mass ratio q = A'/ms/Mwd 
(see Hachisu, Kato, & Nomoto 1999a, for more details). 
Here we modified equation (21) of Hachisu et al. (1999a) 
to include the effect of Roche lobe overflow from L3 point. 
Then the stripping rate is estimated as 



where 



Mstrip = ClMwind, 

VeS ■ 9{q) f v'^a 

Ci = 



')L3-(t>MS \2GM 



(12) 
(13) 



The efficiency of conversion is assumed to be ryoff = 1. 
Here, we assume further that Mv^td = 1-25 Mq, Mms = 
3.0 Mq. Then, we have the separation of a = 4.375i?0, the 
total mass of the binary M — 4.25Mq, and the difference 
of Roche potentials (/)l3 — ^ms = 0.3, and g{q) = 0.03 for 
the mass ratio q = Mms/Mwd ^ 3 (Hachisu et al. 1999a). 
Substituting these values to equation (13), we have 

ci « 0.1 ( ~ 5, (14) 

V600 kms-V 



for the wind velocity of Wwind = 4000 km s^^. Here, the 
maximum wind velocity observed is 4000 km s~^ in the op- 
tical bright state (e.g., Koch et al. 1986; Wood & Lockley 
2000), which is a similar wind velocity in RX J0513. 

5.3. Mass transfer attenuated by winds 

In such a wind phase, the net mass accretion rate to the 
WD, Mace, is modified as 



M^ 



( or e), 



MS 



strip* 



for Mms < Mstrip 
for Mms > Mstrip 



(15) 



Here we adopt a small value of e = 1 x 10 "^Mq yr ^ when 
Mms < Mstrip because the mass accretion to the WD does 
not stop abruptly but probably continues at least for a 
draining time of the accretion disk. We do not know the 
exact draining time of the accretion disk after the mass 
transfer from the MS stops. Alternatively, we just assume 
a small rate of the mass accretion e to mimic the draining 
of the disk during that Mms < Mgtrip- 

To know when the rapid mass accretion resumes, we 
monitor the level of flow as 

^Mflow = Mms - Mstrip - Mace, (16) 

with the initial value of Mflow = (see Fig. 3 of Hachisu & 
Kato 2003b, for more details). This corresponds to when 
the surface of the MS companion fills its critical Roche 
lobe. We regard that the mass transfer is going on when 
Mflow > (overflow level). During that Mstrip > Mms, 
Mflow decreases to a large negative value. Here, we re- 
gard that the mass transfer stops when Mflow < 0. Then 
the mass transfer from the MS stops or its rate drops to 
e. Once the mass accretion rate drops below the criti- 
cal rate, i.e.. Mace < Md-, the wind gradually weakens 
and eventually stops. Then, the stripping effect vanishes 
(Mstrip = 0) and the level of flow (Mflow < 0) begins to 
rise toward Mflow = 0. We start again the rapid mass ac- 
cretion when the level reaches Affiow = 0. This period is 
considered to be a thermal recovering time of the envelope 
of the MS companion. It should be noted that the level 
remains to be Mflow = if a steady-state is reached, i.e., 

Mms - Mstrip - Mace = 0. 

5.4. Time evolution of winds 

Time evolution of the photospheric radius and temper- 
ature of the WD are calculated from Kato & Hachisu's 
(1994) optically thick wind solutions. It is written, after 
Hachisu & Kato (2001b), as 

^ AMenv = Mace " M„uc " Mwind , (17) 

where AMonv is the envelope mass. Mace the mass accre- 
tion rate, Mnuc the decreasing rate by nuclear burning, 
and Mwind the wind mass loss rate of the WD. The photo- 
spheric radius, temperature, and luminosity are given as 
a function of the envelope mass. When the photospheric 
radius shrinks smaller than the binary size, the irradiation 
effects both of the accretion disk and of the companion 
star become important. We include these irradiation ef- 
fects together with a shadowing effect by the disk on the 
companion star as already described in the preceding sec- 
tion (see, e.g., Hachisu & Kato 2001b; Hachisu, Kato, & 
Schaefer 2003; Hachisu & Kato 2003b). 
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5.5. Delay of mass transfer via accretion disk 

An important thing is that equation (17) stands on the 
WD but equations (11) — (16) are on the MS. Therefore, 
the mass accretion rate, Mace, may not be the same on 
each side at the specific time, t, because it takes a viscous 
timescale of the accretion disk for gas to reach the WD 
surface from the MS surface. Then, we have 



Macc(t) 



WD 



Ms^ccit tvis) 



MS 



(18) 



Here we adopt a viscous timescale of 

'^disk -^disk -Rdisk 

V 

the a parameter of Shakura 



^vis — 



40 



/OvisV 

V 0.1 / 



days, (19) 



where ly is the viscosity, oivii 
& Sunyaev (1973), Cs the sound speed, and H the vertical 
height of the accretion disk. We adopt iidisk = 1-4 Rq 
and H/Rdisk ~ 0.1 for Tdisk ~ 30000 - 50000 K and 
Cs ~ 30 km s~^. 

5.6. Model parameters 

Our long-term light curves of V Sge are calculated for a 
given set of parameters: (1) the white dwarf mass, Mwd, 
(2) the companion mass, Mms, and its surface temper- 
ature, TMs,org, (3) the mass transfer rate from the MS 
companion, Mmsi (4) the viscous timescale of the accre- 
tion disk, tyis, (5) the coefficient of stripping effect, ci, 
and (6) the disk parameters of a and (3 and other factors 
which determine the disk shape. The parameters (1) — (5) 
are summarized in Table 1 except for Tms org = 10, 000, 
12,000, and 13,500 K for Mms = 2.5, 3.0,' and 3.5 Mg, 
respectively. As for the disk shapes, we simply adopt the 
same disk shape parameters as shown in Figure 3. Typical 
long-term light curves are shown in Figures 8-11 for var- 
ious white dwarf masses, i.e., Mwd = 1-3, 1.25, 1.2, and 
1.1 Mq, respectively. 

6. NUMERICAL RESULTS OF LIMIT CYCLE MODEL 

6.1. Template light curves 

Now we show our long-term light curves that exhibit 

the transition between the optical high and low states. 
The first case is a binary system consisting of a 1.3 Mq 
WD and a 3.5 Mq MS companion. Assuming Mms = 
20 X IQ-^Mq yr-\ ci = 7.0, and ivis = 41.6 days (i.e., 
81 binary periods), we have followed the time evolution 
of our binary system with the initial envelope mass of 
AMenv = 3 X 10"'' Mq. Figure 8 depicts the obtained 
V^- light curve and color (second panel), accretion rate 
(Afacc on the WD) and wind mass loss rate in units of 
IO^^Mq yr^^ (third panel), radius {Rq) and temperature 
(eV) of the WD photosphere (bottom panel) together with 
the AAVSO observation (top panel). 

At the beginning of time-evolution calculations, the 
mass of the WD envelope is small. As the WD ac- 
cretes matter and its envelope mass increases because of 
Mms = Mg,cc > M^uc- When the envelope mass reaches 
the critical value, AMonv = 4.1 x IQ-'^Mq, the WD pho- 
tosphere expands to i?ph = 0.07 Rq and begins to blow 
a wind. Then the disk surface is blown in the wind and 
the disk size quickly expands from a = 1.15 to a = 6.0 
in dynamical timescales. Here we adopt seven days as the 



transition time, because the massive winds easily blow the 
surface layer of the disk in the wind. The wind mass loss 
rate is about Mwind = 5 x 10~^Mq yr"-'^ only one day af- 
ter the wind starts, which is large enough to completely 
obscure soft X-rays (see, e.g., discussion of Southwell et al. 
1996). The occurrence of optically thick winds can natu- 
rally explain the very rapid emergence and decay of the 
supersoft X-rays accompanied by the transition between 
the optical high and low states (Greiner & van Teeseling 
1998). 

The envelope mass continues to increase and the wind 
mass loss rate also grows. The WD photosphere contin- 
ues to expand and this causes a gradual increase in the 
optical light curve. The mass transfer from the MS is sup- 
pressed when the wind mass loss rate reaches M^ind = 
2.9 x IO-^Mq yr-i, because M^trip = 7Mwi„d = Mms = 
20.0 X 10~^Mq yr~^. However, the mass accretion to the 
WD still continues at a high rate and the wind mass loss 
rate reaches Mwind = 15.0 x 1O~^M0 yr~^ as shown in 
Figure 8 because the mass in the accretion disk is drained 
to the WD for 41.6 days of tvis time delay. The mass ac- 
cretion to the WD eventually stops and the mass of the 
WD envelope begins to decrease due to wind mass loss and 
nuclear burning. The WD photosphere gradually shrinks 
from i?wD,ph = 2.2 Rq to 0.07 Rq during the wind phase. 
This causes the gradual decrease in the optical light curve 
until the end of high state. When the mass of the WD 
envelope decreases to AMenv = 4.1 x 10"'' M©, the wind 
stops. The duration of the wind phase is about 170 days 
as shown in Figure 8 and summarized in Table 1. 

The WD photosphere begins to shrink quickly after the 
wind stops. It takes 12 days to collapse from 0.07 Rq to 
0.03 Rq. The pliotospheric temperature increases quickly 
from 30 eV to 40 eV during this period. It takes 27 days 
to further contract from 0.03 Rq to 0.015 Rq. The photo- 
spheric temperature further increases from 40 eV to 60 eV 
during this period. When the wind stops, the wind mass 
loss rate decreases from M^ind = 1-5 x 10~^Mq yr^^ to 
zero in a day. This means a sharp emergence of supersoft 
X-rays in a day or so as observed in the LMC supersoft 
X-ray source RX J0513.9-6951 (Reinsch et al. 2000). 

During the massive wind phase, the MS companion is 
stripped off by the wind and its surface becomes smaller 
than the critical Roche lobe, i.e., the level of mass flow, 
Mflow, in equation (16) goes to much below the level of 
zero. When the wind stops, the surface of the MS com- 
panion is still smaller than the Roche lobe. It gradually 
recovers up and eventually reaches the level of zero. It 
takes ^ 90 days in Figure 8. It takes a further 41.6 days 
until the mass accretion to the WD resumes, i.e., ivis delay 
time. It takes more 3 days that the WD envelope expands 
to blow a wind again. The total duration of low states 
amounts to 90 -f 40 -|- 3 « 130 days. Thus, the system 
repeats the cycle mentioned above, i.e., a 170 days optical 
high state and a 130 days optical low state. 

The modeled system reaches a limit cycle after one cy- 
cle of high and low states. Thus, the long-term light curve 
modulations are roughly reproduced. The time-averaged 
mass transfer rate to the WD is 3.0 x 1O~^M0 yr~^ 
and 80% of the transferred matter is blown in the wind, 
so that the MS loses its mass at the average rate of 
17.0 x 1O~^M0 yr~^ by stripping effects. These val- 
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ues satisfy the conditions of Mms = -^acc + -^strip and 
-^strip = 7 Mvjind in averaged values. 

Robertson et al. (1997) divided the long-term opti- 
cal brightness of V Sge into three levels, i.e.. the high 
{V < 11.4), intermediate (11.4 < < 12.0), and low 
(y > 12.0) states. In this three-level classification of the 
brightness, we attribute the phase of largely expanded WD 
photospheres (i?wD,ph ^ 0.5 Rq) to the real optical high 
{V < 11.4) state. When the WD photosphere expands as 
largely as > 1 Rq, it contributes more to the UV/optical 
light compared with the disk irradiation. The temperature 
of the WD photosphere is as high as ~ 100, 000 K, while 
the disk has lower temperatures of 10,000 — 50,000 K. 
Thus, the B — V color goes blucward as shown in Figures 
8-11. 

Lockley, Eyres, & Wood (1997) detected radio emission 
from V Sge and estimated the wind mass loss rate of 

M^ind - 8 X 10-« f ^V^' { ..ZT" -i l 
V3 kpcy V1500 km s V 

(20) 

The wind mass loss rate in our model (e.g., Fig. 8) reaches 
-^wind ~ 15 X 1O~^M0 yr~-^ and its averaged value during 
the wind phase is as large as M^ind ^ 5 x 1O~^M0 yr~^, 
which is very consistent with the radio observation. 

Thus, we arc able to reproduce the basic observational 
features (1) — (3) of V Sge summarized in §1. Furthermore, 
our model naturally explains the other observational fea- 
tures: (4) the WD photosphere continues to expand during 
the rapid mass accretion phase and still goes on until ~ 40 
days after the wind started. This makes a visual magni- 
tude peak of ~ 0.5 — 1.0 at the head of optical high states. 
Large expansions (~ 2 Rq) of the WD photosphere are 
very consistent with the photospheric radius (~ 1.7 Rq) 
of the WD in the optical bright state shown in Figure 5. 

6.2. White dwarf mass 

We have examined other WD masses, i.e., Mwd = 1-25, 
1.2, and 1.1 Mq, with the companion mass being kept to 
be Mms = 3.5 Mq. The results are shown in Figures 9, 
10 and 11. In general, more massive white dwarfs have 
a shorter duration of the wind phase. The durations of 
optical high and low states are summarized in Table 1 for 
various white dwarf masses and other parameters. 

Our guiding principle is that (1) we tune ci and Mms 
to satisfy the duration of the low state (about 130 days), 
because wc do not know the exact number of ci, but it 
may be around 3—10 from equation (14), and that (2) we 
also tune Mms to reproduce the optical peak shape in the 
bright state. Each tuned parameters are summarized in 
Table 1. 

For the case of 1.25 Mq WD + 3.5 Mq MS, we choose 
the parameters of ci = 8.0 and Mms = 25 x 10~^Mq yr~^. 
Then, the durations of the wind phase (optical high state) 
and of the no wind phase (optical low state) are about 220 
and 130 days, respectively. The observational light curve 
shape of long-term variations is well reproduced by this 
1.25 Mq WD + 3.5 Mq MS model. 

When the rapid mass transfer from the MS companion 
stops, we still keep a low rate of mass accretion to the WD, 
i.e., e in equation (15). To check the effect of e, we double 
the value of e and calculate the long-term light curve. The 



duration of the optical low state becomes a bit shorter, 
120 days, as listed in Table 1. On the other hand, the 
duration of the optical high state becomes a bit longer, 
230 days, with the total period being almost not changed, 
350 days. Wc may conclude that the effect of e is not so 
large as long as its value is small compared with the rapid 
mass transfer rate. This also examines the effect of the at- 
mospheric mass transfer, although the real number of the 
mass transfer rate coming from the tenuous atmosphere of 
the MS companion is not estimated, but it probably does 
not much exceed 1 x 10"'' Mq yr~^ (e.g., Ritter, Zhang, & 
Kolb 2000). 

For the case of 1.2 Mq WD + 3.5 Mq MS, we have the 
parameters of ci = 10.0 and Mms = 30 x IO^^Mq yr-^ 
Then, the durations of the wind phase (optical high state) 
and of the no wind phase (optical low state) are about 
290 and 140 days, respectively. The calculated duration 
of the optical high state (wind phase) is too long to be 
compatible with the observation. 

For the case of 1.1 Mq WD 3.5 Mq MS, we choose the 
parameters as ci = 14.0 and Mms = 40 x IO^'^Mq yr~-^. 
The durations of the optical high state and of the optical 
low state are about 430 and 150 days, respectively. In this 
case, too, the calculated duration of the optical high state 
is much longer than that of the observation. 

When the durations of the low states arc set to be similar 
to those for Mwd = 1.3 Mq, the duration of high states 
becomes much longer in the less massive white dwarfs. 
From the theoretical light curve shapes and the durations 
of optical high states, we may conclude that the WD mass 
in V Sge is somewhere between 1.3 and 1.25 Mq. In what 
follows, we use Mwd = 1-25 Mq as already adopted in 
equation (1). 

6.3. Mass transfer rate of MS 

The mass transfer rate of the MS companion depends 
on its evolutionary stage and binary parameters (e.g., Li 
& van den Heuvel 1997; Langer et al. 2000). Since we 
do not know the exact number of the mass transfer rate 
in V Sge, we have examined various values of the original 
mass transfer rate Mms, i-e., 0.64, 1.0, 2.0, 5.0, 10.0, 20.0, 
25.0 (Fig. 9), and 40.0 in units of IO'^Mq yr"! for the 
1.25 Mq WD + 3.5 Mq MS. Their results are summarized 
in Table 1 and the long-term light curves are shown in Fig- 
ure 12. Here, we fix the other parameters, i.e., fvis = 41.6 
days and ci = 8.0. The duration of optical high states 
becomes longer for the higher mass transfer rates of Mms • 
On the other hand, the duration of the optical low state 
hardly changes. 

As the mass transfer rate of Mms increases, an optical 
peak at the head of the high state becomes more promi- 
nent. This is because the WD photosphere expands to as 
large as 1 — 2 Rq and both the irradiation effect by the 
disk and the WD brightness itself increase and contribute 
much to the y- light. From this peak luminosity, we adopt 
the mass transfer rate of Mms = 25 x 10~^Mq yr~^ for 
the case of 1.25 Mq WD + 3.5 Mq MS. 

6.4. Effect of viscous timescale 

The viscous timescale of accretion disks is not theoret- 
ically determined yet. Therefore, we examine the depen- 
dency of the value of t^is on our modeled light curves. 
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Here, we adopt four other different values of tvis = 13.9 
days (27 rotation periods), fvis = 27.8 days (54 rota- 
tion periods), tvis = 55.5 days (108 rotation periods), and 
tyis = 69.4 days (135 rotation periods), as summarized in 
Table 1 and shown in Figure 13. The other parameters 
are fixed, i.e., Mms = 25 x IQ-'^Mq yr"^ and a = 8.0 for 
the binary model of 1.25 Mq WD + 3.5 Mq MS. 

For the case of tvis = 13.9 days, we have no low states 
because the rapid mass accretion resumes 13.9 days after 
the condition Mms > .^strip is satisfied. This time is be- 
fore the wind stops. For the case of tvis = 27.8 days, we 
have a short duration of low states, 30 days, much shorter 
than the observational ones. 

On the other hand, for a relatively long timescale of 
tvis = 55.5 days, we have a rather long duration of low 
states, 240 days, much longer than the observational ones. 
This is because the rapid mass accretion still continues 
during 55.5 days after the condition Mms < Mgtrip is ful- 
filled and, as a result, a massive and fast wind strips off 
the surface of the MS companion deeply, and it takes a 
much longer time to recover the level of zero in equation 
(16). For the case of tvis = 69.4 days, we have a very long 
duration of low states, 340 days. The steady hydrogen 
shell-burning stops at 270 days in the low state, that is, 
70 days before the optical high state starts. 

We have a reasonable duration of optical low states when 
the viscous timescale is around tvis = 40 days, which is 
roughly corresponding to ctvis = 0.1 of Shakura & Sun- 
yaev (1973) as seen in equation (19). 

6.5. Mass stripping effect 

It is a difficult work to accurately estimate the effect of 
mass-stripping, that is, to determine the coefficient of ci 
in equation (12). Therefore, we calculate other four cases 
of ci, i.e., ci = 3.0, ci = 5.0, ci = 10.0, and ci = 12.0 and 
examine how the light curves depend on the coefficient. It 
is obvious that strong modulation of the mass transfer rate 
never occurs for the case of ci < 1 because there is no case 

of Mms < Mstrip = ciMwi„d(< Mwind < M^cc < Mms)- 
Our calculated models and their parameters are listed in 
Table 1. 

For the case of ci = 3.0 shown in Figure 14. we have 
no optical low states, mainly because the mass-stripping 
effect is small and, as a result, the recovery time of the MS 
envelope is rather short. For the case of ci = 5.0, the du- 
ration of optical low states is about 10 days, still because 
of small mass-stripping effect. 

On the other hand, for the case of ci = 10.0, the mass- 
stripping effect increases and the recovery time of the MS 
envelope becomes drastically longer. If we further increase 
the value of ci = 12.0, the duration of low states reaches 
280 days, where the steady hydrogen shell-burning stops 
at 270 days in the low state, i.e., 10 days before the optical 
high state starts. 

6.6. Flat segment to active segment 

In the long-term light curve behavior, there are peri- 
ods with/without strong optical modulations. Simon & 
Mattel (1999) called them active/flat segments, respec- 
tively. The active segments can be understood as a limit 
cycle of the optical high and low states as already shown 
in the previous subsections. Here, we interpret the nature 



of the flat segments. Simon et al. (2001) pointed out 
that the flat segments finished by another decrease into a 
real low state and further suggested that the intermediate 
(medium) state is classified as the flat segment (see §6.1). 

Our WD wind model provides a natural explanation of 
such behavior if the mass transfer rate from the MS com- 
panion varies in time together with the mass-stripping ef- 
fect. In the calculation shown in Figure 15, we gradually 
increase Mms from 1 x lO'^M© yr'^ to 25 x IO-^^Mq yr'^ 
in 1200 days and, at the same time, strengthen the mass- 
stripping effect from ci = 1.0 to 8.0. This may imply that 
the wind velocity (the mass-stripping effect) is relatively 
slow (small) at the low wind mass-loss rate regime but 
gradually increases proportionally to the wind mass-loss 
rate. No modulations appear when the mass stripping ef- 
fect is as small as ~ 1.0 — 1.5 and the mass transfer rate 
is also as low as ~ 1 x 10~^Mo yr~^. This figure shows 
that modulation appears in the light curve as both the 
mass transfer rate and the mass stripping effect increase 
in time. 

6.7. Dependence on the secondary mass and the 
inclination angle 

The long-term light curve behavior depends also on the 
other various parameters. In this subsection, we exam- 
ine the long-term light curves by changing the other sys- 
tem parameters. (1) The light curve also depends on the 
mass of the MS companion. We have calculated two other 
masses of Mms = 2.5 M© and 3.0 M© with the inclina- 
tion angle oi i = 77°. The results are shown in Figure 
16. The overall light curve features are not changed but 
the brightness becomes fainter in the less massive compan- 
ion masses. (2) We have also checked the dependency of 
the light curve on the inclination angle, i.e., i = 82° and 
i = 60° as shown in Figure 16. The depth of the optical 
low state becomes much deeper at the lower inclination 
angles, that is, much deeper at i = 60° than at i = 80°. 

7. DISCUSSION 

7.1. The He II 4686 emission line in the high state 

Phase-resolved spectra in the optical high state have 
been observed by several groups (Gies, Shaffer, & Wiggs 
1998; Lockley & Wood 1998; Robertson, Honeycutt, & 
Pier 1997; Wood & Lockley 1997). Lockley (1998) pre- 
sented high-quality, phase-resolved spectra of He II A4686 
and H/3 emission lines in the optical high state (see also, 
Lockley & Wood 1998; Patterson et al. 1998; Wood & 
Lockley 1997). He showed six kinematic components in 
the He II A4686 emission line. Gies et al. (1998) summa- 
rized them as follows: 

(1) A broad, flat-topped component with steep edges that 
extend further to the red than the blue. This broad fea- 
ture also dominates Ha, and double-Gaussian flttings of 
the extreme wind radial velocity variations are similar in 
both He II and Ha (Gies et al. 1998; Robertson et al. 
1997). The phase-resolved behaviors of these two sharp 
edges, which are attached by number "1", are shown in 
Figure 17. The left dashed-line corresponds to the blue 
edge while the right dashed-line does to the red edge. 

(2) A central, narrow, nebular line. This component ap- 
pears stationary in both He II and Ha. The Ha profiles 
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appear doublc-pcakcd in both Gies et al. (1998) and 
Robertson et al. (1997), because of the presence of both 
He II A6560 and Ha. The Ha component indicates an out- 
flow with a velocity of ~ 200 km s~^. The phase-resolved 
behavior of this Unc is also shown in Figure 17, to which 
number "2" is attached. 

(3) An outer blue peak that accelerates blueward after 
primary eclipse and attains minimum velocity near the 
extreme blue wind just prior to phase 0.5, to which num- 
ber "3'' is attached. This feature is also scon in Ha. 

(4) An outer red peak, to which number "4" is attached, 
best seen near phase 0.7, that moves along with the red 
sharp edge of the broad component 1. 

(5) An inner red peak that appears near phase 0.6, to 
which number "5" is attached, moves redward from line 
center, and then disappears again near phase 0.1. This 
feature is seen in Ha. 

(6) An inner blue peak that emerges from the central neb- 
ular line to follow the motion of star 1 (primary compo- 
nent), to which number "6" is attached. This is seen as 
a slight blueward extension of the Ha peak in the radial 
velocity range to -1-200 km s~^ between phase 0.2 and 
0.5. 

With the above features, we presume three line forma- 
tion regions of the He II A4686 in our WD optically thick 
wind model: the first one is the WD massive wind, the 
second one is the irradiated surface of the disk, and the 
third one is the irradiated side of the MS companion and 
the matter which is stripped off from the MS companion. 
The first and second regions are denoted by filled arrows 
and short outlined arrows, respectively, in Figure 18. 

The red sharp edge (component 1 and 4) almost fol- 
lows the motion of star 1 (the primary WD) with ~ 
250 km s~^ but red-shifted by ~ -1-1500 km s~^. It sug- 
gests that this component is originated from the WD wind 
with the velocity of ^ 1500 km s"-'^. The sharp blue edge 
of emission component 1 is not consistent with this WD 
wind model because it does not show a sinusoidal feature 
of Ki ~ 250 — 300 km s~^. This feature may be explained 
by the effect of self-absorption by the WD wind itself, al- 
though it is heavily contaminated by the C III+N III com- 
plex in the blue edge over —1000 km s~^. The spectra 
of Lockley & Wood (1998) and Patterson et al. (1998) 
certainly show a weak feature of the C III-I-N III complex. 

Emission components 3 and 5 also have a similar Ki ~ 
250 — 300 km s^^ and follow the motion of star 1 with a 
smaller red shift of ~ -1-600 km s~^ for component 5 and 
a blue shift of ~ —600 km s~^ for component 3. This 
suggests the velocity of the disk surface flow as fast as 
~ 600 km s~^ (see also Fig. 4). 

For Ha emission, Gies et al. (1998) obtained a ra- 
dial velocity curve of Ki ~ 380 km s~^ by double- 
Gaussian fitting for the inner flat-topped structure with 
a = 1000 km s~^. This implies a bit faster disk surface 
outflow, which is reasonably explained if the line forma- 
tion region of Ha is further out than that of He II 4686 
and the disk surface flow there is more accelerated. 

The central emission components 2 and 6 (dash-dotted 
lines in Fig. 17) show anti-phase motions of the other com- 
ponents. This feature indicates that components 2 and 6 
are originating from star 2. As we have already discussed, 
the MS companion (star 2) loses matter stripped by the 



WD wind in the optical high state. This gas outflow is 
much slower compared with the WD wind and the disk 
surface flow, may be as slow as ~ 100 — 200 km s~^. 

7.2. The He II 1640 emission line in the low state 

The main features of the observed He II 1640 spectra 
in the optical low state (Wood & Lockley 2000) arc sum- 
marized as follows: (1) At phase zero the line is roughly 
symmetrical about its center although an excess in the red 
peak is seen. (2) As the binary phase increases, the emis- 
sion on the blue side of the line increases until at phase 
~ 0.3 with emission out to a velocity of ~ —1500 km s~^. 
At this phase a sharper red edge appears at ^ 300 km s~^. 
(3) At phase 0.4 the line is roughly symmetric again. (4) 
Then, from phase 0.5 to 0.8, the situation is opposite to 
phase from 0.1 to 0.4, with the extended emission now 
occurring on the red side out to ^ 1600 km s~^ and the 
sharper blue edge at ^ —200 km s^^. 

Wood & Lockley (2000) explained the above features 
by the colliding winds. However, these features can be 
qualitatively explained even by our WD model. Fukuc & 
Hachiya (1999) showed that a disk surface flow is driven 
up to 3000 — 4000 km s~^ when the central WD is as lumi- 
nous as the Eddington luminosity. The disk surface flow 
is line-driven, therefore, we see an optically thin wind go- 
ing outward along the disk surface. When the disk edge 
is asymmetric as illustrated in Figure 19, the disk surface 
flow is also asymmetric because the elevated edge of the 
disk suppresses effective acceleration. (1) At phase zero, 
we can see both the elevated disk edge and the extended 
emission region in the angle between phase 0.1 and phase 
0.4 (see Fig. 19). The emission region of the elevated disk 
edge contributes to an excess at the red side of the central 
peak. (2) As the binary phase increases, the blue side of 
the line becomes more prominent. During phase 0.1 to 0.4, 
the disk motion is toward us and its velocity is increasing 
to —300 km s~^ at phase 0.25. This orbital motion makes 
the sharp red edge of the line blueward as shown in Figure 
4 of Wood & Lockley (2000). (3) At phase 0.4 the line 
is roughly symmetric but with an small excess in the blue 
side. This is very consistent with the observed one. (4) 
From phase 0.5 to 0.8 the behavior of the line is opposite 
to that of phase 0.1 to 0.4 as easily understood from Figure 
19. Since the disk goes away from us during this phase, 
the sharp blue edge of the line moves redward. This is also 
consistent with Figure 4 of Wood & Lockley (2000). At 
phase 0.9 the line is almost symmetric again but with a 
very small excess in the red side and this is also consistent 
with the observed line profiles. 

7.3. The C IV doublet 1548 and 1551 in the low state 

The observational features of the C IV spectra in the 
optical low states (Wood & Lockley 2000) are summa- 
rized as follows: (1) At phase 0.5 — 0.7 the blue absorption 
feature belonging to the blue member of the C IV dou- 
blet becomes largest with the P Cygni absorption reaching 
~ —700 km s~^. (2) As the phase increases the extent de- 
creases, reaching only ~ —250 km just prior to eclipse. 
(3) The extent of the absorption feature increases again af- 
ter eclipse to phase 0.5 but the absorption itself is rather 
weak. 

To explain the above behavior we assume that the line 
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formation region of C IV doublet is much more extended 
compared with that of He II 1640. An illustrated configu- 
ration of our model is shown in Figure 20. In the direction 
of phase 0.4 to 0.9, the disk surface flow goes up along 
the elevated edge of the disk involving the disk matter 
dragged. Therefore, the disk surface flow crosses the line 
of sight and absorbs the C IV lines. On the other hand, 
in the direction of phase 0.9 to 0.4 the disk surface flow 
virtually does not cross the Hue of sight (see Fig. 20). 
As a result, the absorption is very weak. The disk sur- 
face flow has a large angle against the line of sight, for 
example, 30 — 60° or so. Then, the extension of the the 
P Cygni absorption feature is limited at a smaller velocity 
like ~ —700 km s~^ even if its terminal velocity is as large 
as ~ 1000- 1500 km 8"^ 

7.4. Other emission features in the low state 

Phase-resolved spectra of the He II A4686 and H/3 in 
the optical low states were obtained by Diaz (1999). The 
radial velocities of the emission lines were measured by 
a double-Gaussian convolution method. Their results are 
summarized as follows: 

(1) The radial velocity amplitudes of He II A4686 and H/3 
agree with each other, i.e., K ~ 230 km s^^. by Gaussian 
components half-separated by 800 km s~^ for He II A4686 
and by 1100 km s'^ for H/3. 

(2) Positivc-to-negative velocity crossing occurs at phase 
0.93, that is, the largest positive velocity is reached at 
phase 0.68 and the smallest negative velocity is attained 
at phase 0.18. 

(3) The H/? line flux presents a sharp maximum during 
phase 0.0, which indicates that the Balmcr line forma- 
tion region is not eclipsed, while the He II A4686 line flux 
reaches minimum at phase 0.0, indicating a partial eclipse 
of a highly ionized region. 

We attribute the He II A4686 line formation region to 
a relatively inner portion of the disk but the H/3 line for- 
mation region both to an outer portion of the disk and 
the irradiated MS companion (see Fig. 19). The disk 
surface flow is accelerated by the central WD and its out- 
flow velocity gradually increases outward. This tendency 
is roughly consistent with the Gaussian component veloc- 
ities of 800 km s~^ (inner portion of the He II A4686 line) 
and 1100 km s~^ (outer portion of the H/3). 

The disk structure is asymmetric as shown in Figure 19, 
so that we observe the largest positive velocity at phase 
0.65 (midpoint between phase 0.4 and 0.9) and the low- 
est negative velocity at phase 0.15 if the disk surface flow 
is much faster than the orbital motion. Adding the or- 
bital motion of the WD to the disk outflow, the phase of 
the largest positive velocity should be shifted a bit later 
because the orbital motion reaches its largest velocity at 
phase 0.75. Therefore, phase 0.68 is reasonable. For the 
lowest negative velocity, the phase should also be delayed 
from phase 0.15 because the orbital motion of the WD 
component reaches the lowest negative value at phase 0.25. 
In this sense, phase 0.18 is also reasonable. 

During phase 0.5 to 0.8, the blue shifted component of 
the disk flow is blocked by the elevated disk edge as seen 
in Figure 19, so that an extended blue wing is probably 
cut in the spectra as well as the self-absorption feature. 
Diaz (1999) showed contour maps of phase-resolved line 



intensities both for the He II A4686 and H/3 lines. The 
above feature is clearly shown in Figures 4 and 5 of Diaz 
(1999) both for the He II A4686 and H/3 hues. 

Closely looking at the radial velocity change of the 
reddest contour in Figures 4 and 5 of Diaz (1999), 
we easily find that the radial velocity reaches maximum 
of ~ 1500 km at phase ^ 0.7 and minimum of 
~ 900 km s~^ at phase ~ 0.1, which suggests an am- 
plitude of Ki ~ 300 km s~^ and an outflow velocity of 
^ 1200 km s^^. These kinematic velocities are consistent 
with our WD wind model. However, the bluest contour 
shows the minimum velocity of —1000 km s~^ at phase 0.2 
and the maximum velocity of —500 km s~^ at phase 0.6 
for the He II A4686 line. This asymmetric feature is very 
similar to those of He II A4686 emission line in the optical 
high state (see Fig. 17), although the outflow velocity is 
smaller in the low state. It is therefore likely that the blue 
edges of these emission lines arc self-absorbed even in the 
optical low state, unless it is signiflcantly contaminated by 
C III-I-N III complex in the blue edge over —1000 km s~^. 

To summarize, our WD wind model possibly gives us 
natural explanations of the observational features, i.e., 
He II A4686 emission (together with Ha) in the optical 
high states and He II A1640 emissions, C IV doublet A1548 
and 1551, and He II A4686 emission (together with H/3) in 
the optical low states, at least qualitatively. 

7.5. Weak soft X-ray flux 

Diaz & Steiner (1995) pointed out that soft X-ray flux 
of V Sge is too weak to be compatible with the typical 
supersoft X-ray sources, for example, at least two or three 
orders of magnitude lower than that of CAL 87. In this 
subsection, we try to resolve this problem. Our numeri- 
cal results are summarized as follows: (1) the distance to 
V Sge is 3 kpc, (2) the temperature of the WD is between 
30 — 60 eV as seen in Figure 9, and (3) the bolometric 
luminosity is 1.7 x 10'^^ erg s^-'^. If wc adopt the color 
excess of E{B — V) = 0.3 and the corresponding hydro- 
gen column density of A^h ~ 2 x 10^^ cm~^ (Gorenstein 
1975), the ROSAT/PSPC count rate is calculated to be 
larger than 2 counts s^^. However, the observed rate is 
only 0.02 counts (Greiner & van Teeseling 1998), at 
least, two orders of magnitude smaller than that of our 
WD model. 

The central WD is always occulted by the elevated disk 
edge in the CAL 87 model (model d) of Schandl et al. 
(1997). In this model, a substantial part of soft X-rays 
can penetrate the elevated edge of the disk, i.e., the spray 
is semi-transparent at short wavelength (sec private com- 
munication of F. Meyer in Lederle & Kimeswenger 2003). 
Very recently, Suleimanov, Meyer, & Meyer-Hofmeister 
(2003) formulated X-ray radiation reprocessing in the el- 
evated edge of the disk. In their results, some fraction of 
soft X-rays leak from the elevated edge that consists of 
many small size cloud blobs and hot (^ 5 x 10^ K) in- 
tercloud medium. The mass transfer rate in CAL 87 is 
estimated to be as small as ~ 1 x 1O~^M0 yr~^ and two 
orders of magnitude smaller than that for our V Sge model. 
It is very likely that the density of spray is much (e.g., two 
orders of magnitude) larger in V Sge than in CAL 87. This 
much more opaque elevated disk edge results in a very low 
flux of soft X-rays even if the system configuration (orbital 



V Sagittae 



11 



period and inclination) is almost the same between these 
two systems. 

7.6. Galactic counterpart of RX J0513.9- 6951 

Greiner & van Teeseling (1998) concluded that the 

model suggested for RX J0513. 9—6951 cannot explain the 
X-ray and optical variabilities of V Sgc. Their reasons 
are as follows: (1) both the transitions, from the faint to 
bright and the bright to faint, are very rapid, so that they 
may occur in within one day (compared to the smooth 
dcchnc of several days in RX J0513.9-6951). (2) The ex- 
pected optical eclipse would become deeper in the brighter 
state, opposite to what has been observed (Patterson et al. 
1998). 

Their assumed model for RX J0513.9— 6951 is based on 

the photospheric expansion/contraction model of the WD 
envelope (e.g., PakuU et al. 1993; Reinsch et al. 1996, 
2000). Their conclusions suggest only that the photo- 
spheric expansion/contraction model does not apply to V 
Sge. As already shown in Hachisu fc Kato (2003a, b), the 
photospheric expansion/contraction model cannot explain 
the X-ray variability even for RX J0513.9— 6951 because 
the emergence/decay timescale of supersoft X-rays is also 
as short as ^ 1 — 2 days (Reinsch et al. 2000). 

On the other hand, our model can naturally explain the 
above two observational features as follows: (1) the irra- 
diated disk is a main source of the optical light and its 
flaring-up edge is very variable in a dynamical timescale, 
i.e., a day or so. This also naturally explain the optical 
variability from night to night (Lockley et al. 1999; Mader 
& Shafter 1997). The difference in the optical variabilities 
between V Sge and RX J0513.9— 6951 may come from the 
difference in the inclination angle, i.e., z = 70 — 80° for V 
Sge but i = 20-30° for RX J0513.9-6951 (Hutchings et al. 
2002; Hachisu & Kato 2003b). The irradiated area of the 
flaring-up disk viewed from the Earth varies much more 
largely in the binary system with a higher inclination an- 
gle. (2) As recently reported in Cowley et al. (2002), the 
eclipse depth of RX .J0513.9— 6951 is deeper in the optical 
low state than in the optical high state. This is the same as 
what has been observed in V Sge. Our models can repro- 
duce the orbital modulations both in the optical high and 
low states both for V Sge and RX J0513.9-6951. Thus, 
we may conclude that V Sge is essentially the same type of 
objects as RX J0513.9— 6951, and its Galactic counterpart. 

7.7. Relevance to Type la supemovae 

Thus, our optically thick wind model of mass- accreting 
white dwarfs can reproduce many observational features 
of V Sge. We suggest that V Sge is the second example 
of accretion wind evolution, which is a key evolutionary 
process to Type la supernovae in a recently developed evo- 
lutionary scenario (Hachisu et al. 1996, 1999a,b; Hachisu 
& Kato 2001b). In the accretion wind evolution phase, 
the donor transfers mass to the white dwarf at a high rate 
of ~ 1 X 1Q~^Mq yr~^ or more. In an old picture, the 
white dwarf expands to a giant size to form a common en- 
velope and results in merging. However, the white dwarf 
begins to blow a massive wind in such a high mass ac- 
cretion rate instead of simply expanding. A formation of 
common envelope is avoided and the separation of the bi- 
nary is almost unchanged. In the accretion wind evolution 



phase, the white dwarf accretes matter and burns hydro- 
gen atop the white dwarf core. Then, the white dwarf can 
grow in mass at the critical rate, Mcr 1 x 10~^Mq yr"-'^, 
up to the Chandrasekhar mass limit and explodes as a 
Type la supernova. With the first and second examples 
of accretion wind evolution RX J0513.9— 6951 (Hachisu & 
Kato 2003a, b) in the LMC and V Sge in our Galaxy, we 
may conclude that accretion wind evolution is established 
in the evolutionary scenario to Type la supernovae. In 
other words, accretion wind evolution commonly occurs in 
the supersoft X-ray sources when the mass transfer rate 
exceeds the critical rate of ~ 1 x 1O~^M0 yr~^. 

The life time of V Sge phenomena is as short as ^ 10'' yr, 
because the mass transfer rate is as large as Mms ~ 
1O~^M0 yr~^. V Sge belongs to the helium-rich supersoft 
X-ray source channel in the Type la supernova scenario of 
Hachisu et al. (1999b). This channel can produce about 
one; Type la supernova per millennium, indicating that 
we have a chance to observe about a hundred V Sagittae 
stars in our Galaxy. Steiner & Diaz (1998) have already 
listed four V Sagittae stars in our Galaxy and discussed 
their similar properties. Although the mass of the com- 
panion star to the WD is not clearly identified yet, their 
orbital periods fall in the range of 0.2 — 0.5 days, being very 
consistent with the orbital periods that the new scenario 
predicted (see Fig. 3 of Hachisu & Kato 2001b). 

8. CONCLUSIONS 

We have proposed a model of long-term optical vari- 
ations in the peculiar binary V Sge based on an opti- 
cally thick wind model of mass- accreting white dwarfs. 
Our binary model, which is consisting of a mass-accreting 
white dwarf (WD), an accretion disk around the WD, and 
a lobe-filling main-sequence (MS) companion star, repro- 
duces orbital light curves both for the optical high and low 
states. Moreover, our WD wind model naturally explains 
the transition between the optical high and low states. 

Here, we summarize our main results: 

1. When the mass accretion rate to the white dwarf ex- 
ceeds a critical rate, the white dwarf photosphere expands 
to blow a fast and massive wind. The disk surface is 
blown in the wind and optically thick parts of the disk 
extend to a large size over the Roche lobe. The largely 
extended disk, which is irradiated by a WD photosphere, 
contributes much to the optical light. This corresponds to 
the optical high states. 

2. The fast and massive wind hits the donor star and 
strips off the surface layer. The mass transfer rate from 
the donor is attenuated. As a result, the wind gradually 
weakens and stops. Then the disk shrinks to a size of 
the Roche lobe. This (no wind) phase corresponds to the 
optical low states because the disk irradiation area also 
becomes small. 

3. In the optical bright state, we are able to reproduce 

the orbital light curves of V Sge with a large extension of 
the disk. In the optical faint state, the orbital light curves 
are reproduced with a small size (similar size to the Roche 
lobe) of the disk having an elevated edge as seen in the 
supersoft X-ray sources. 

4. The massive wind easily shields supersoft X-rays. 
Therefore, soft X-rays are not detected in the wind phase, 
i.e., in the optical high state. This is consistent with the 
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soft X-ray behavior found by Grcincr & van Tccscling 
(1998). Furthermore, our wind model also explain the 
hard X-ray component during the optical high state if 
they are shock-origin (Greiner & van Teeseling 1998). 

5. The observed radio flux indicates a wind mass loss rate 
of ~ IO-^Mq yr-i (Lockley et al. 1997). Our optically 
thick wind model of mass-accreting white dwarfs naturally 
explain such a high mass loss rate. 

6. Various features of the emission lines both in the op- 
tical high and low states are also naturally explained by 
our model. 

7. Optically thick winds from mass-accreting WDs play 

an essential role in a recently developed evolutionary 
scenario of Type la supernovae (Hachisu et al. 1996, 
1999a, b). This process is called accretion wind evolution. 
V Sge is the second example of accretion wind evolution 



with the first example of the LMC supersoft X-ray source 
RX J0513.9-6951 (Hachisu & Kato 2003a,b). With these 
two examples of accretion wind evolution, we may con- 
clude that optical high/low states with wind being on/off 
are a common physical feature in the supersoft X-ray 
sources if the mass transfer rate exceeds the critical rate 
of ~ 1 X IO-^Mq yr-i. 

8. V Sge is a strong candidate for Type la supernova 
progenitors. 

We thank V. Simon for providing us machine readable 
data of the AAVSO observation. This research has been 
supported in part by the Grant-in-Aid for Scientific Re- 
search (11640226) of the Japan Society for the Promotion 
of Science. 
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^mass accretion rate Afacc is modulated by winds 

'^mass accretion rate is slightly attenuated but not modulated largely by winds, so that it is almost constant 
in time 
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MS companion accretion disk 

Fig. 1. — A schematic configuration of our V Sge system in the faint state is illustrated. The main-sequence (MS, left) companion is 
strongly irradiated by the shell- burning white dwarf (WD, right) so that the irradiated hemisphere mainly contributes to the light from the 
MS companion. If the stronger fluorescent O III line comes from the surface of the MS companion, it certainly traces closely the irradiated 
hemisphere rather than the non-irradiated side. As a result, the orbital velocity observed is lower than that of the MS companion itself. 
Crosses around the MS surface denote emission of the stronger fluorescent O III line and C.G. indicates the center of gravity of the binary 
system. The azimuthal angle ip is measured from the x-axis of orbital phase 0.5 (the secondary eclipse phase). 
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Fig. 2. — The orbital velocities of each component are plotted against the MS companion mass for three inclination angles, i.e., i = 80°, 70°, 
and 60° (from upper to lower), when the WD mass is MwD = 1.25 Mq. If we adopt Ki = JfwD ~ 300 — 340 km s~^, the mass of the MS 
companion ranges from M2 = Mms = 3.0 Mq to 3.5 Mq and the velocity of the MS companion ranges from K2 = Kms = 130 km to 
110 km s-l. 
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a) optically thick massive wind phase 




b) static (no nnossive wind) phase 



Fig. 3. — Configurations of our V Sge model are illustrated: (a) in the massive wind phase (optical high/soft X-ray off), and (b) in the 
static shell-burning phase, i.e., the rapid mass accretion phase after the wind stops (optical low/soft X-ray on). The cool component (right) is 
a main-sequence (MS) companion (S.OMq) filling up its inner critical Roche lobe. The north and south polar areas of the cool component are 
irradiated by the hot component (1.25 Mq white dwarf, left). The separation is a = 4.375ii0; the effective radii of the inner critical Roche 
lobes are R* = 1.34i?0, and R2 = R2 = 2.0Rq, for the primary WD and the secondary MS companion, respectively. The surface layer of the 
disk is blown in the wind, like a free stream. This disk surfa<;e stream is optically thick near the original disk region but becomes optically 
thin far outside the original disk due to geometrical dilution effect. We regard the transition place from optically thick to thin as the outer 
edge of the disk in the massive wind phase (see also Fig. 4a). This edge extends over the MS companion as shown in the figure. The disk 
shrinks to a normal size (here 1.15 times the Roche lobe size) in several orbital periods after the wind stops. A rapid mass accretion resumes 
to make a spray around the disk edge in more several orbital periods. The inclination angle is i = 77° and the aaimuthal angle is rp = 80°. 
Here, the azimuthal angle is measured from the angle at the secondary eclipse, i.e., at the phase with the WD component in front of the MS 
companion. 
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WD photosphere 



(a) optically thick massive wind phase 




(b) no massive wind (static) phase 



Fic;. 4. — Optically thick parts of the disk arc illustrated: (a) in the massive wind phase (optical high/soft X-ray off), and (b) in the static 
shell-burning phase (optical low/soft X-ray on). See text for more details. 
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Fig. 5. — Calculated V light curves and B — V colors are plotted against the binary phase (binary phase is repeated twice from —1.0 to 
1.0) together with the observational points (Mader & Shafter 1997; Herbig et al. 1965). Solid lines denote our calculated models. Various 
system parameters are shown in the figure together with the best fitted parameters, (a) the bright (optical high) state, (b) the faint (optical 
low) state. 
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(b) static (no massive wind) phase 



Fig. 6. — Same as Fig. 3 but for the eclipse minimum (0 = 180°); (a) in the massive wind phase (optical high/soft X-ray ofT) and (b) in the 
static shell-burning (no massive wind) phase after the wind stops (optical low/soft X-ray on). The WD surface is completely blocked by the 
MS companion but a part of the disk is seen from the Earth. 
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(a) optically thick massive wind phase 




(b) static (no nnossive wind) 



phase 



Fig. 7. — Same as Fig.3 but for the secondary minimum (</> = 0°): (a) in the massive wind phase (optical high/soft X-ray off) and (b) in 
the static shell-burning (no massive wind) phase after the wind stops (optical low/soft X-ray on). The WD surface is blocked by the edge of 
the disk. 
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Fig. 8. — Top panel: (a) AAVSO visual data are plotted against time (JD 2,448,000+). Second panel: (b) V-magnitude light curves (thick 
solid) of our V Sge model (1.3 Mq WD and 3.5 Mq MS) are plotted together with the calculated color {B — V)c (thin solid). The calculated 
y-light curve is made by connecting the brightness at orbital phase 0.4. The inclination angle is assumed to be i = 77°. Third panel: (c) 
Mass accretion rate to the WD (Maco thick solid) and wind mass loss rate from the WD (M^inj, thin solid), both in units of 10~^Mq yr~^. 
Bottom panel: (d) Photospheric radius of the WD envelope in units of Rq (thick solid) and surface temperature of the WD envelope in units 
of eV (thin solid). The model parameters are summarized in Table 1. 
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Fig. 10.— Same as Figure 8, but for 1.2 Mq WD + 3.5 Mg MS. 
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Flo. 12. — Calculated l/-magnitudc light curves (thick solid) of V Sge for various mass accretion rates are plotted against time (JD 
2,448,000+), together with the calculated color (B — V)c (thin solid). The mass accretion rates are sliovim in each panel. The other 
parameters are fixed and summarized in Table 1. Here, the binary is 1.25 Mq WD + 3.5 Mq MS. The inclination angle is assumed to be 
i = 77°. (g) AAVSO data of V Sge tciken from Simon & Mattei (1999) are also plotted. 
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Fig. 13. — Same as Fig.l2, but for various values of the delay by the accretion disk, tvis- The value of t^is is shown in each panel. 
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Fig. 14. — Same as Fig.l2, but for various values of the stripping effect, ci. The value of ci is shown in each panel. 
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Flc;. 15. — Transition from a fiat segment to an active segment in tfie optical light curve is calculated. Both ci and Mms a-re gradually 
increasing with time. The value of ci is going up from ci = 1.0 to ci = 8.0 and the mass transfer rate is also growing from Mms = 
1 X lO-'^M© yr-i to Mms = 25 x lO-'^M© yr'^, both in 1200 days. 
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Fig. 16. Same as Fig. 12, but for other two cases of the companion mass and inchnation angle, i.e., (a) Mms = 2.5 Mq with the inclination 
angle of i = 77°, (b) Mms = 3.0 Mq and i = 77°, (c) AAVSO visual data taken from Simon & Mattel (1999), (d) Mms = 3.0 Mq and 
i = 82°, and (e) Mms = 3.0 Mq and i = 60°. The shapes of V- light curves are not so different from each other, although the depth of the 
low state becomes deeper for the lower inclination angle (i = 60° ) . 
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Fig. 17. — Velocities for He II 4686 emission components are illustrated against the orbital phase. Attached numbers denote the He II 4686 
emission components defined by Lockley (1998). These velocities are compiled from Fig. 1 of Lockley & Wood (1998) and Figs. 1 and 2 of 
Wood & Lockley (1997). See text for more details. 
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Fig. 18. — Line formation regions of He II A4686 in the optical high sates are schematically illustrated. The WD (left component) blows 
a massive wind, the velocity of which is as fast as 2000 — 2500 km s^-"- (long filled arrows). The surface layer of the disk is dragged outward 
and its surface extends up to several times the Roche lobe size. Then, the disk surface is accelerated to 1000 km (short outlined arrows). 
Motions of binary components are denoted by open arrows, which axe starting from the center of each component. The binary orbital phase 
is also denoted by outer open arrows. 
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Fig. 19. — An optically thin, weak disk surface flow (wind) is driven by the central hot white dwarf when the white dwarf is very luminous 
(Pukue k, Hachiya 1999). (a) The optically thick photo-surfaces (photospheres) are shown in a bird view, (b) The formation region of He II 
1640.476 emission line is illustrated (enclosed by a thick solid line). Arrows in the line formation region indicate the velocities of the disk 
surface flow. The disk surface flows can easily flow out from phase 0.9 to 0.4 and are possibly accelerated up to 1500 km s~^ or more (Pukue 
& Hachiya 1999). On the other hand, it cannot be accelerated so much between phase 0.4 to 0.9 because the elevated edge of the disk 
hampers the acceleration. The outer arrows represent the orbital phase (the primary eclipse is phase 0) viewed from the Earth. 
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Fig. 20. A cross sectional view of the disk surface flow (wind) is shown. No or weak self-absorption (P Cygni profile) of the C IV doublet 
is observed from the orbital phase of 0.1 while a rather strong self-absorption is seen from the orbital phase of 0.6. We assume that the line 
formation region of the C IV doublet is more extended compared with that of the He II 1640 emission in Figure 19. 



